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Abstract 

The reaction pp — > Pf(r]rj)p s has been studied at 450 GeV/c. For the first time a par- 
tial wave analysis of the centrally produced 7777 system has been performed. Signals for the 
/o(1500), /o(1710) and /2(2150) are observed and the decay branching fractions of these states 
are determined. 
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The 7777 , 777/ and rfr}' channels are important for the determination of the gluonic content of 
mesons, since it is thought likely that glueballs will decay with the emission of 77s and 77's [p], |2[. 
In addition, central production is proposed as a good place to produce glueballs via Double 
Pomeron Exchange (DPE) || |3|. Hence the interest of the present paper which studies the 
centrally produced 7777 system. 

The r/77 channel has been studied previously in radiative J/ip decays by the Crystal Ball 
experiment ||, in n~p interactions by the NA12 experiment 0, H|, in central production by 
NA12/2 [0] and in pp annihilations by the Crystal Barrel Collaboration || and E760 at Fermi- 
lab |§. In all these reactions evidence for the /q(1500) has emerged. In addition, in pp annihila- 
tion evidence is claimed for a / (2100) @, [Kj] and an 7777 decay mode of the / 2 (1950) [nj. In n~p 
interactions the NA12 experiment claim a /o(2000) ||. In radiative J/ip decays evidence was 
reported for the /j(1710) H. In central production the NA12/2 experiment claimed evidence 
for a / 2 (2175) 

In this paper a study is presented of the 7777 final state formed in the reaction 

PP-*Pf(w)Ps (1) 

at 450 GeV/c. It represents more than a factor of 6 increase in statistics over the only other 
data on the centrally produced 1777 final state |7[ and moreover will present the first partial wave 
analysis of this channel in central production. The data come from the WA102 experiment which 
has been performed using the CERN Omega Spectrometer, the layout of which is described in 
ref. | TT|| . Reaction ([!]) has been isolated using the following decay modes: 



77 — > 77 77 — > 7T + 7r _ 7r 
77 — > 7T + 7r _ 7r 77 — > 7T + 7r~7r° 

The above decay modes account for 38.8 % of the total. 

Fig. |I|a) shows a scatter plot of M(77) versus M (77) which has been extracted from the 
sample of events having two outgoing charged tracks and four 7s reconstructed in the GAMS- 
4000 calorimeter using momentum and energy balance. A clear signal of the 7777 channel can be 
observed. Fig. [lfb) shows the 77 mass spectrum if the other 77 pair is compatible with being 
an 77 (0.48 < M(77) < 0.62 GeV) where a clear 77 signal can be observed. The 7777 final state 
has been selected using the mass cuts described above. 

Fig. [l|c) shows a scatter plot of M(77) versus M(7r + 7r~7r°) for the sample of events having 
four outgoing charged tracks and four 7s reconstructed in the GAMS-4000 calorimeter after 
imposing momentum and energy balance. A clear signal of the 7777 channel can be observed. 
Fig. |l|d) shows the 7r + 7r~7r° mass spectrum if the 77 mass is compatible with being an 77 
(0.48 < M(77) < 0.62 GeV) and fig |l|e) shows the 77 mass spectrum if the 7r + 7r~7r° mass is 
compatible with being an 77 (0.52 < M(7r + 7r~7r°) < 0.58 GeV). The 7777 final state has been 
selected by requiring that 0.48 < M(77) < 0.62 GeV and 0.52 < M(7r+7r-7r°) < 0.58 GeV. 

Fig. |TJf ) shows a scatter plot of M(tt + 7t~7t°) versus M(7r + 7r~7r°) for the the sample of events 
having six outgoing charged tracks and four 7s reconstructed in the GAMS-4000 calorimeter 
after imposing momentum and energy balance. A signal of the 7777 channel can be observed. 
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Fig. shows the tt + 7t~7t° mass spectrum if the other 7r + 7r~7r° combination has a mass com- 
patible with being an rj (0.52 < M(7r + 7r~7r°) < 0.58 GeV). The 7777 final state has been selected 
using the mass cuts described above. 

The background below the rj signal has several sources including combinatorics, fake gam- 
mas and other channels. The combinatorial background is removed, in part, in the selection 
procedure. The remaining background varies from 16 % to 32 % dependent on the decay topol- 
ogy. Three methods have been used to determine the effects of this background; studying the 
side bands around the rj signal, studying events that do not balance momentum and using 
event mixing. All three methods give a similar background which peaks near threshold in the 
rjrj mass spectrum. To illustrate this point, superimposed on figs. |l]b), d), e) and g) as a shaded 
histogram are the respective mass distributions when the central system has a mass greater 
than 1.3 GeV. As can be seen the background below the rj signal is reduced with respect to 
the total sample. In the remainder of this paper the method used to determine the background 
will be the one using events that do not balance momentum. The resulting distributions well 
represent the background in the data. 

The 7777 mass spectra from each decay mode are very similar and the combined mass spec- 
trum is shown in fig. 0a) and consists of 3351 events. Superimposed on the mass spectrum as 
a shaded histogram is the estimate of the background. Fig. |b) shows the background sub- 
tracted mass spectrum. The mass spectrum has a threshold enhancement and has peaks at 
1.5 and 2.1 GeV and a shoulder in the 1.7 GeV region. Superimposed on fig ||b) as a shaded 
histogram is the expected contribution from the /2(1270) (47 events) and ^2(1525) (3 events) 
using the observed ^(1270) signal in the 7r + 7r~ final state |12|], the observed ^2(1525) signal 
in the K + K~ final state [T3| and correcting for the experimental acceptance and using the 
PDG [[14] branching ratios. As can be seen the majority of the signal at 1.5 GeV is not due to 



the $(1525). 

A Partial Wave Analysis (PWA) of the centrally produced rjrj system has been performed 
assuming the rjrj system is produced by the collision of two particles (referred to as exchanged 
particles) emitted by the scattered protons. The z axis is defined by the momentum vector of 
the exchanged particle with the greatest four-momentum transferred in the rjrj centre of mass. 
The y axis is defined by the cross product of the momentum vectors of the two exchanged 
particles in the pp centre of mass. The two variables needed to specify the decay process were 
taken as the polar and azimuthal angles (6, <fi) of one of the 77s in the rjrj centre of mass relative 
to the coordinate system described above. 



The acceptance corrected moments v47rt iM , defined by 

J(fi) = ^t L0 V L °(fi) + 2 £ t LM Re{Yl ! (n)} 

L L,M>0 

have been rescaled to the total number of observed events and are shown in fig. [| The moments 
with M > 2 and all the moments with L > 4 (e.g. the £44 and t^o moments shown in fig. ||) are 
small and hence only partial waves with spin I = and 2 and absolute values of spin z-projection 
m < 1 have been included in the PWA. 



The amplitudes used for the PWA are defined in the reflectivity basis |T5| . In this basis the 
angular distribution is given by a sum of two non-interfering terms corresponding to negative 
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and positive values of reflectivity. The waves used were of the form with J = S and D, 
m = 0, 1 and reflectivity e = ±1. The expressions relating the moments an d the waves 

( J^J are given in ref.JIIJ. Since the overall phase for each reflectivity is indeterminate, one wave 



in each reflectivity can be set to be real (Sq and for example) and hence two phases can 
be set to zero (<p s - and 4> D + have been chosen). This results in 6 parameters to be determined 
from the fit to the angular distributions. 

The PWA has been performed independently in 80 MeV intervals of the f]7] mass spectrum. 
In each mass an event-by-event maximum likelihood method has been used. The function 

N 

F = -J2 ln{I(n)} + t LM e L M (2) 

i=l L,M 

has been minimised, where N is the number of events in a given mass bin, e^u are the effi- 
ciency corrections calculated in the centre of the bin and tLM are t ne moments of the angular 
distribution. The moments calculated from the partial amplitudes are shown superimposed as 
histograms on the experimental moments in fig |3|. As can be seen the results of the fit reproduce 
well the experimental moments. 

The system of equations which express the moments via the partial wave amplitudes is 
non-linear which leads to inherent ambiguities. For a system with S and D waves there are 
two solutions for each mass bin. In each mass bin one of these solutions is found from the fit 
to the experimental angular distributions, the other one can then be calculated by the method 
described in ref. |T5[. In the case under study the bootstrapping procedure is trivial because 
the Barrelet function has only two roots (Z± with i = 1,2) and their real and imaginary parts 
do not cross zero as functions of mass, as seen in fig ||]c) and d). In order to link the solutions 
in adjacent mass bins, the real parts of the roots are sorted in each mass bin in such a way that 
the real part of the first root should be larger than the real part of the second root (real parts 
of the two roots have different signs) . For the first solution the imaginary parts of both roots 
are taken positive, the second solution is obtained by complex conjugation of one of the roots. 

In this case two different PWA solutions are found. One solution is dominated by the 
S-wave the other solution has the events split between the different D-waves. By definition 
both solutions give identical moments and identical values of the likelihood. A Monte Carlo 
study has shown that if the input distribution is really dominated by a D-wave (i.e. Dq , 
or Di) then both solutions of the PWA will be dominated by that D-wave. However, if the 
input distribution is dominated by the S-wave then two different solutions will result. One of 
the PWA solutions will be dominated by the S-wave, the other solution will be split equally 
between the D-waves and hence the physical solution can still be determined. The physical 
solution is shown in fig. |j. The S'J'-wave for the physical solution is characterised by a broad 
enhancement at threshold and a peak at 1.5 GeV. A broad enhancement is also seen in the 
.D^-wave at 2.1 GeV. Superimposed on the waves as a histogram is the result of running the 
PWA on the background events. 

Fig. |5| shows the background subtracted Sq and Dq waves. The PWA analysis has also been 
performed by extending the likelihood function given in equation (U) to include the background 
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subtraction, namely 



N N bg 

F = -J2 ln{I{ty} + ]T t LM e LM + ]T ln{I(n)} 

i=l L,M i=l 

where iV& s is the number of background events. Since the background is dominantly S-wave the 
results of this method are similar to those shown in fig. |5|. 

The S^-wave has been fitted using a K-matrix parameterisation similar to that used to 



fit the K + K~ spectrum [17] with the addition of an incoherent background term. Poles 
have been introduced to describe the / (1370), /o(1500) and /o(1710) with parameters fixed 
to those found from the coupled channel fit to the 7r + 7r~ and K + K~ spectra [|l7j], namely 
M(/ (1370)) = 1312 MeV, r(/ (1370)) = 218 MeV, M(/ o (1500)) = 1502 MeV, r(/ (1500)) = 
98 MeV and M(/„(1710)) = 1727 MeV, r(/ (1710)) = 126 MeV. The fit describes well the data 
but, due to the size of the errors, it is not possible to conclude about the need for a further 
resonance above 2 GeV. The inclusion of the /o(1710) is essential to describe the spectrum. 

The Dq wave has been fitted using two spin 2 relativistic Breit-Wigners and a linear back- 
ground. The first Breit-Wigner is used to describe the /2(1270) with mass and width fixed 



to the PDG values [jl4| and the second to describe the peak at 2.1 GeV where we have used 
M = 2130 MeV, T = 270 MeV which are the parameters found for / 2 (2150) observed in the 
K + K~ final state ||13|| . As can be seen the fit well describes the data. This state is compatible 
with the /2(2175) previously observed by the NA12/2 experiment in the r]f] final state J7|. 

The error bars introduced by the partial wave analysis do not allow the parameters of the 
resonances to be determined from a free fit to the waves. Instead, we have performed a fit to 
the total mass spectrum shown in fig. |2|b) using an incoherent sum of the expressions used to 
fit the Sq and Dq waves. The parameters for the /o(1370) and /2(1270) have been fixed to 
the values used above. The fit is shown superimposed on fig. |2|b) and for the scalar resonances 
yields sheet II T-Matrix poles fl8|] at 



/o(1500) M = (1510 ± 8) -i(55± 8) MeV 

/o(1710) M = (1698 ± 18) -i ( 60 ± 13) MeV 



These parameters are consistent with the PDG WM values for these resonances and our previous 



fits 0. For the / 2 (2150) the fit gives 

/ 2 (2150) M = 2151 ± 16 MeV, T = 280 ± 70 MeV. 



The parameters found are compatible with those from the K + K channel [13 



After correcting for the unseen decay modes and using data from the previously observed 
mi [O, [16] final state the branching ratio r/rj /tttc has been calculated for the /o(1500) from the 



fit to the S wave and gives : 

/o(1500) -> rjrj 
/o(1500) -> 7T7T 



0.18 ±0.03 



This value agrees well with the value that can be derived from the PDG 1T3] of 0.23 ± 0.10. 
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For the / 2 (1270) the ratio is 

/ 2 (1270) -> rrn 
/ 2 (1270) -> 7T7T 

which is compatible with the PDG |L4| value. 



(3±1) x 10" 



For the /q(1500), /q(1710) and /2(2150) the branching ratio r\r\ /KK has been calculated 
from the WA102 data to be: 



/o(1500) - 


-> T)T) 


/o(1500) - 


> KK 


/o(1710) - 


-> ?7?7 


/o(1710) - 




/ 2 (2150) - 


-> ?7?7 



0.48 ±0.15 
0.78 ±0.14 



/ 2 (2150) -> KiT 

The branching ratios for the /o(1500), /o(1710) and / 2 (2150) differ from that of a known ss 
state, the /2(1525), which has a ratio of 0.12 ± 0.04 which is consistent with the SU(3) 
prediction for such a state [jnj. No signal is seen for the /2(1525) in the rjr] final state of WA102 
and an upper limit for its decay to this channel has been calculated, which gives 

£(1525) -VV =<QU 



£(1525) -> KK 
which is compatible with the PDG [Of value given above. 



For the /o(1370) there is considerable uncertainty due to the background subtraction which 
mainly affects the Sq wave at threshold. The strongest / (1370) signal has been observed in 
the 47r final state therefore the branching ratio rji] /47r has been calculated to be 

/o (1370)^ 

/o(1370) -> 4vr 1 ; 

which is compatible with the Crystal Barrel measurement of (1.7 ± 0.9) xl0~ 3 [^TJ . 



In the D-waves there is no evidence for the / 2 (1950) which has been claimed to have been 
seen in the rjr] final states formed in pp annihilations ||10|| . Since the / 2 (1950) is clearly seen 
in the 47r final state of experiment WA102 [f2TJ, an upper limit for its decay to rjrj has been 



calculated and gives 

/ 2 (1950)- W<50xio _ 3 (go% CL) 



/ 2 (1950) -> 4vr 

Hence this would imply that if the observation in the 7777 final state of pp annihilations is 
correct, then a very large signal for the / 2 (1950) should be seen in the 47r final state of the 
same experiment. 



The /o(1500) has previously been observed in the rnr JT2|, lg], KK |]13| , An [^1| and r/rj' [p2 
final states of the WA102 experiment. The relative decay rates for the /o(1500) are calculated 
to be: 

7T7T : KK : rjr) : rjrf : 4tt = 1 : 0.33 ± 0.07 : 0.18 ±0.03 : 0.096 ± 0.026 : 1.36 ±0.15 
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The /o(1370) is below 777/ threshold. The remaining relative decay rates for the /o(1370) are: 

7T7T : KK : 7777 : 4tt = 1 : 0.46 ± 0.19 : 0.16 ± 0.07 : 34.0+f 

No signal was observed for the /o(1710) i n the 47r |^T[ or the 777/ channels Therefore, an 
upper limit has been calculated for its decay to these final states. For the /o(1710) the relative 
decay rates are: 

7T7T : KK : 7777 : 1717' : 4tt = 1 : 5.0 ± 0.7 : 2.4 ± 0.6 :< 0.18 (90 % CL) :< 5.4 (90 % CL) 

In summary, a partial wave analysis of the 7777 channel has been performed for the first time 
in central production. Clear evidence is found for an 7777 decay mode of the / (1500), / (1710) 
and /2(2150) and the decay branching ratios of these states have been determined. 
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Figures 



Figure 1: Selection of the 7777 final state. For the reaction pp — > PfPsi^l) a) M(77) versus 
M(77) and b) M(77) if the other 77 pair is in the 77 band (0.48 < M(77) < 0.62 GeV). For 
the reaction pp — > pfp s (ir + 7i^7i 2'-f) c) M(77) versus M(7r + 7T~7r°), d) M(7r + 7r~7r°) if the 77 
mass is in the 77 band (0.48 < M(77) < 0.62 GeV) and e) M(77) if the 7r + 7r~7r° mass is in 
the 7] band (0.52 < M^Tr^Tr ) < 0.58 GeV). For the reaction pp -> p f p s (n + n-n + n-n°n°) f) 
M(7T + 7r _ 7r°) versus M(7r + 7r~7r°) and g) M(7r + 7r~7r°) if the other it + it~it mass is in the 77 band 
(0.52 < M(7T + 7r~7r°) < 0.58 GeV). Superimposed as a shaded histogram are the respective 
histograms when the central system has a mass greater then 1.3 GeV. 



Figure 2: a) The 7777 mass spectrum. Superimposed as a shaded histogram is an estimation 
of the background contribution, b) The background subtracted 7777 mass spectrum with fit 
described in the text. Superimposed as a shaded histogram is an estimation of the ^(1270) 
and /2(1525) contribution. The c) Real and d) Imaginary parts of the roots (Z, see text) as a 
function of mass obtained from the PWA of the 7777 system. 



Figure 3: The \f^ixt LM moments from the data. Superimposed as a solid histogram are the 
resulting moments calculated from the PWA of the 7777 final state. 



Figure 4: The physical solution from the PWA of the 7777 final state. Superimposed as a shaded 
histogram is an estimation of the background contribution. 



Figure 5: The background subtracted a) Sq and b) Dq waves with fits described in the text. 
The parameters used to describe the resonances are those previously determine from the nn 
and KK channels. 
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Figure 2 
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Figure 5 
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